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New summation inequalities and their 
applications to discrete-time delay systems 

Le Van Hien and Hieu Trinh 


Abstract 

This paper provides new summation inequalities in both single and double forms to be used in stability 
analysis of discrete-time systems with time-varying delays. The potential capability of the newly derived 
inequalities is demonstrated by establishing less conservative stability conditions for a class of linear 
discrete-time systems with an interval time-varying delay in the framework of linear matrix inequalities. 
The effectiveness and least conservativeness of the derived stability conditions are shown by academic 
and practical examples. 


Index Terms 

Summation estimates, discrete-time systems, time-varying delay, linear matrix inequalities. 

I. Introduction 

It is well known that time-delay frequently occurs in practical systems and usually is a source of 
bad performance, oscillations or instability lUl, 121. Therefore, the problem of stability analysis and 
applications to control of time-delay systems are essential and of great importance for both theoretical 
and practical reasons which have attracted considerable attention, see, for example ISl- ITOl and the 
references therein. 

Among existing works which concern with stability of linear time-delay systems, the Lyapunov- 
Krasovskii functional (LKF) method plays an essential role in deriving efficient stability conditions. 
Based on a priori construction of a Lyapunov-Krasovskii functional combining with some bounding 
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DRAFT 


techniques 0|, Q, ifTOl - lfT^ . improved delay-dependent stability conditions for continuous/discrete time- 
delay systems were derived in terms of tractable linear matrix inequalities Q, lITOl . ifT^ . llT3l . However, 
the design of such Lyapunov-Krasovskii functional and especially the techniques used in bounding the 
derivative or difference of constructed Lyapunov-Krasovskii functional usually produce an undesirable 
conservatism in stability conditions. Therefore, aiming at reducing conservativeness of stability conditions, 
it is relevant and important to improve some fundamental inequalities to be used in establishing such 
stability criteria ifTOl . |[T4l . 

Note that, most of the aforementioned works have been devoted to continuous-time systems. Besides, 
the problem of stability analysis and control of discrete-time systems with time-varying delay is very 
relevant and therefore it should be receiving a greater focus due to the following practical reasons. 
Firstly, with the rapid development of computer-based computational techniques, discrete-time systems 
are more suitable for computer simulation, experiment and computation. Secondly, many practical systems 
are in the form of nonlinear and/or non-autonomous continuous-time systems with time-varying delays. 
A discretization from continuous-time systems leads to discrete-time systems described by difference 
equations which inherit the similar dynamical behavior of the continuous ones ITSl . In addition, the 
investigation of stability and control of discrete-time systems requires specific and quite different tools 
from the continuous ones. Thus, stability analysis and control of discrete-time delay systems have received 
considerable attention in recent years |[T6l - l[2^ . Most recently, novel summation inequalities were derived 
ll25l . ll2^ by extending the Wirtinger-based integral inequality Q. These summation inequalities provide 
a powerful tool to derive less conservative stability conditions for discrete-time systems with interval 
time-varying delay in the framework of tractable linear matrix inequalities. 

In this paper, new summation inequalities which provide an efficient tool for stability analysis of 
discrete-time systems with time-varying delay are fisrt derived. Inspired by the approaches proposed in 
IHOI . |[T4l for the continuous-time systems, new summation inequalities in both single and double forms 
are derived by refining fhe discrefe Jensen inequalifies. If is worfh nofing fhaf fhe obfained resulfs in fhis 
paper fheoretically encompass fhe summation inequalities proposed in 1251 . ||2^ . Furthermore, unlike 1251 . 
1261 . we prove that the proposed inequalities do not depend on the choice of first-order approximation 
sequences. By employing these new inequalities, a suitable Lyapunov-Krasovskii functional is constructed 
and less conservative stability conditions are derived for a class of discrete-time systems with interval time- 
varying delay. To illustrate the effectiveness of the proposed stability conditions, an academic example 
and a practical satellite control system are provided. These examples show that our stability conditions 
provide significant improvement over existing works in the literature. 
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The paper is organized as follows. Section II presents some preliminary results. New summation 
inequalities and their applications to stability analysis of a class of discrete-time systems with interval 
time-varying delay are presented in Section III and Section IV, respectively. Numerical examples to 
demonstrate the effectiveness of the obtained results are also given in Section IV. 


II. Preliminaries 

Notations: Throughout this paper, we denote Z and Z+ the set of integers and positive integers, 
respectively, M” the n-dimensional Euclidean space with vector norm ||.||, the set of n x m real 


matrices. For matrices A,B€ 


A 0 


\ co\{A,B} and diagj^d, i3} denote the block matrix 


and 


, respectively. A matrix P G is positive (negative) definite, write P > 0 (P < 0) if 

0 PJ 

x^Px > 0 {x'^Px < 0) for all x G M", x 7 ^ 0. We let denote the set of symmetric positive definite 
matrices. For any A G R”'^"', He(A) stands for A + . For a,b ^ a < b, Z[a, b] denotes the set of 
integers between a and b. For a sequence u : Z[a, 6 ] M”', we write Uk = u{k), k G Z[a, b], A denotes 
the forward difference operator, that means Auk = Uk+i—Uk- For any two sequences u, u : Z[a, 6 ] ^ M”, 
it is obvious that UkAvk = A{ukVk) — Ufc+iAufc. 

The following inequalities which are widely used in the literature can be easily derived by using Schur 
complement lemma. 

Lemma 1: (Jensen’s inequalities) For a given matrix R G integers b > a, any sequence u : Z[a, b] 

R”, the following inequalities hold 

'^ulRuk>-(^Ukf R{'^Uk), ( 1 ) 


k=a 


k=a 


k=a 
h k 




i(i + l) 

k=as=a k=as=a 

where i = b — a + 1 denotes the length of interval [a, b] in Z. 


( 2 ) 


k=as=a 


III. New summation inequalities 

In this section, new summation inequalities are derived by refining ([T]), (ID). In the following, let us 
denote Ji{u) and J 2 {u) as the gap of ([T]) and ([ 2 ]), respectively, that is the difference between the left-hand 
side and the right-hand side in ([1]) and ©.By refining ([U and (O, we aim to find new lower bounds 
for Ji{u), J 2 {u) other than zero. 
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Lemmal: For a given matrix R G S^, integers b > a, any sequence u : Z[a, 6 ] M", the following 

inequality holds 


T / \ ^ + 1) I 5(f' + l)(£ + 2)^ A'Tp/- 

^ - £(£ - 1 ) ^ £(£ - 1)(£2 + 11)^2 ^C 2 


(3) 


b b k 

where (^ = vi - j^V 2 , (2 = vi - j^V 2 + and vi = ^k, V 2 = '^ Y1 '^3 = 

k=a k=as=a 

b k s 

EEEu^. 

k=a s=a i=a 

Proof: Note at first that if = 1 then Ci = C 2 = 0 and thus ([3]) holds. Now assume that i> 1. We 
use the idea of hilevel optimization to get ([3]) hy refining ([T]). To fhis, for a sequence u : Z[a, b] M"', 
we define an approximation sequence n : Z[a, 6] —>• as follows 


1 ^ 

Vk = Uk - ^'^Uk + akXi + /3fcX2 


(4) 


k=a 


where and /3fc are two sequences of real numbers and Xi)X 2 £ are constant vectors which will 

he defined lafer. From dH) we have 

b b b 

^ vlRvk = Ji (u) + 2xf afeUfc) + 2x2 -R(^ /^kUk) 


k=a 


k=a 


k=a 


}C^(^k)xjRvi - -n(^h)xlRvi 


k=a 


k=a 


(5) 


+ C^ak)xlRxi + C^Pk)x^RX2 

k=a k=a 

b 

+ 2{'^akl3k)xlRX2- 

k=a 


k-1 

Let Uk = Ui for k > a, Uk = 0 for k = a, then Uk = and, consequently, akUk = A{akUk) — 

i=a 

Uk+iAak- Taking summation from a to 6 gives 

b b b 

^ ^ OfcMfc — ci 6 +i 'y ^ Uk y ^ Aq;^. ( 6 ) 

k=a k=a k=a 

For any first-order sequence ak which can he written as = co(fe — a) + ci, cq / 0, we have 

06-1-1 = CqI + Cl, Aok = Co, 2 ^ “fc = *20- - -h Ci£. 

k=a 


This, in regard to ®, leads to 


o 2 ^ 

“^XiRC^ akUk) - jC^ ak)xlRvi = ao{£ + l)xrf?Ci- 


(7) 


k=a 


k=a 
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Similar to ® we have 


huk = /3b+i Uk+iAPk- 

k=a k=a k=a 

At this time we define the sequence Uk = 'Wfc+i = and thus 

'^k+l^Pk ^i^Pk'^k') iifc+ 1 ^ f^k- 

For convenience, we choose j3k = {k — a)^ — i{k — a) + then Yl\=al^k = 0, fih+i = 
A/3b+i = ^ + 1 and A'^{(3k) = 2. Note also that 

b b fc+l b fc+l b k 

^uk+l=^^u^=Y^ yz “*=yzyz Ui+l = V3. 


Therefore 


k=a k=as=a k=a s=a+l k=as=a 

^ 2 ^ ^ 1 
‘^X2 f^kUk) - Pk)X2 = 3 X 2 

k=a k=a 


( 8 ) 


where ^3 = — l)r’i — 6(£ + l)r;2 + 12r;3. 

On the other hand, from (@1) and note that Ylk=a h = 0, we readily obtain Yll=a '^k = {Ylk=a Q^fc)xi- 
This, together with (l5]l, dT]) and ([8]), leads to 

1 


Ji{v) = Ji{u) + co{e + i)xi’^Ci + 


+ 




k=a 


k=a 


xJRxi 


(9) 


+ (y^/ 3 fc)X 2 -RX2 + 2 (y^ afc^fc)xT^X2. 


k=a 

It can he verified hy some direct computations that 

b , b 


k=a 


2^^k-j[2^ak} =Co- 
k=a 


k=a 

/?2 1^//J2 


,2 _ £(r - i){e^ + 11) 


12 

b 


180 


) ^ ^ (^kl3k 


= -Co- 


i{e - 1) 


k=a k=a 

By injecting those equalities into ® we then obtain 


12 


Jiiv) = Ji{u) + Coii + Ijxf-RCi + 


cl£{f - 1) 


12 


Xi-RXi + J 


( 10 ) 


where J = 1x2 ^Cs - cq- g ai . igo 


hlRX2. 


6 -xi Rx 2 + 

Now, at the first stage we define xi = —Ci^ where A is a scalar, then by Lemma [U Ji(v) > 0, it 

Cq 

follows from (ITOl ) that 

/ /f/'-l'i „\ ^ 

( 11 ) 


Jl (^^) > (^ + 1) ( A - ClRCl - J- 


12 
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The function A — 


12 


attains its maximum at A = and hence xi = then 




coi{i-iy 


from (fTTI) we obtain Ji(u) > RQi — J- In addition, hy injecting xi = ^ pip_i\ Qi into J we then 


obtain 


Ki=T) 

Jii”) > - tit— 


co^y-iy 


£(£- 1 ) 


180 

{i + l){^ + 2) J. 

-5-A2 -h.C2- 


( 12 ) 


As this stage, we define X 2 = ~30C2^ is a scalar, then by some similar lines when dealing with (fTTI) 
we finally obtain (fS]) which completes the proof. ■ 

Remark 1: The proof of Lemma |2] can be shortened by a specific selection of au, for example, au = 

Remark!: Lemma |2] in this paper generalizes the summation inequality derived in Lemma 2 in ll2^ 
and Lemma 3 in ll25l by the following points. Firstly, the inequality provided in Lemma |2] in this paper 
encompasses both the inequalities proposed in Lemma 2 in Il26l and Lemma 3 in ||251 since a positive 
term is added into the right-hand side of dT). Secondly, and most interesting is that, dT) can be derived 
from the approximation dD for any first-order sequence au = cq/c -)- ci, cq 7 ^ 0 whereas some special 
cases of (fH) were used to derive Lemma 2 in ||26l and Lemma 3 in (251 . Thirdly, a unify approach is 
introduced to derive some new lower bounds of summation estimate in both single and double form 
proposed in Lemma |2] and the following lemmas. 

Lemma3: For a given matrix R E S+, integers b > a, any sequence u : Z[a, 6 ] —)• the following 

inequality holds 

(13) 


Mu) > 


£{£^ - 1 ) 

where U = V 2 - 

Proof: Similar to Lemma O when £ = 1, ^4 = 0 and (fT3l ) is trivial. Assume that £ > 1 . By the same 
approach used in deriving ([3]), we now construct the following approximation 

U U 

2 






EE Us + akX 


(14) 


k=as=a 


for a given sequence u : Z[a, b] 


\ Similar to dS]) 


EE vJRvs = J2{u) + 2x^ fl(EE“.-) 


k=a s—a 

4 


k=as—a 
b k 


(15) 


£{£+l) 




k=aS—a 


k=a s=a 


k=as=a 
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For any first-order sequence ak = co{k — a) + ci, cq / 0, by some similar lines in the proof of Lemma 
| 2 ] we have 

^ 4-rn(f^2) ^ 

(16) 


Tr \ T^ \ , cli{i + — 1) rp 4co(£ + 2) j, 

J2[v) = J2{u) --X Rx 4 -^-X RQa- 


36 - - ' 3 

From Lemma [U J 2 {v) > 0, and by choosing x = “^C 4 ^ it follows from (fT^ that 

J2{u) >{£ + 2) Ua - cJrCa 


(17) 


which yields ([T3]) for A = The proof is completed. ■ 

Remarks: The summation inequalities proposed in Lemma |2] and Lemma |3] in this paper give a new 
lower bound for the gap Ji(m) and J 2 {u) of the discrete Jensen’s inequalities, respectively. In other 
words, new refinements of the celebrated Jensen’s inequalities have been derived in this paper. 

Remark 4: The double summation inequality provided in (fT3l) is closely related to the function-based 
double integral inequality proposed in flOl although the proof (fTSl) is based on a simple idea, refining the 
classical discrete Jensen’s inequality. Furthermore, as shown in the proof of Lemma [3l inequality ([T3l) 
can be derived from (fT4b for any first-order sequence a^. 

Remarks: As discussed in 1261, some coefficients in ([3]) and (fT3l) might be difficult to handle, especially 
in applications to discrete-time delay systems. Therefore (O and (fTSl) will be reduced to simpler forms 
as presented in the following corollary. 

Corollary 1: (Refined Jensen-based inequalities) For a given matrix 7? E S^, integers b > a, any 
sequence u : Z[a, 6 ] —)• the following inequalities hold 

b 


^ulRuk > \vlRvi + - 


k=a 
b k 


1 

Cl 

T 

-1 

CO 

0 

_ 1 


Cl 

1 

1 

_1 


1 - 

0 

o\ 


1 

to 

1_ 

r 

-1 

T r 

-1 


-1 






V2 

Ca 


R 0 
0 8R 


V2 

C4 


(18) 


(19) 


k=as—a 

where £, V2Xi^ C 2 and ^4 are defined in ([3]) and ([T3l) . 

Proof: The proof is straight forward from (O, (fT3]) and thus is omitted here. 


IV. Stability of discrete-time systems with time-varying delay 

This section aims to demonstrate the effectiveness of our newly derived summation inequalities through 
applications to stability analysis of discrete-time systems with interval time-varying delay. 
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A. Stability conditions 


( 20 ) 


Consider a linear discrete-time system with interval time-varying delay of the form 

x{k -f 1) = Ax{k) + A(ix{k — h{k)), k > 0, 
x{k) = (t){k), t G [-/i2,0], 
where x{k) G M"" is the state, A, A^ G are given matrices, h{k) is time-varying delay satisfying 

hi < h{k) < /i2, where hi < /i2 are known positive integers. For simplicity, hereafter the delay h{k) 
will he denoted hy h. 

Let {e*}i<j<io he the row basic of and e* = e* (8* In- We denote A = {A — In)ei -f AdC^ and 

x{k) 

I x{k — hi) 


Co{k) = col { 


ai{k) = 

= 

U5{k) = 
ue{k) = 


1 


x{k — h) 

x{k - /l2) 
k 


ui{k) 

k'2{k) 

i-'sik) 


Tihi) 

1 


x{s),U 2 {k) = 


U4{k) 

v^{k) 

VQ{k) 

1 


k—hi 


s=k—hi 

k—h 


Tih2 - h) 

1 

7 (/r - hi) 

1 


l{h2 - h) 
T{h)=h + l,y{h) = 


^ x{s),Vi{k) = 
s=k—h2 

—hi k—hi 

s=—h i=k-\-s 
—h k—h 

Y. Y 

s=—h 2 i=k-\-s 

T{h)T{h + l) 


T{h-hi) 
1 


Y 


lihi) 


s=k—h 
0 k 

Y Y 

s=—hi i=k-\-s 


Q{h) = col{ei,r(/ii)e5,r(/i - hi)ee + T{h2 - h)e7,y{hi)e8} 


fli = col{-7l, 62,63 -f 64, T(/ri)e5} , 

0,2 = Co 1{0,6 i ,62 -f 63,r(/li)6i}, 

Fi = Col{6l - 62,61 -f 62 - 265,61 - 62 “f 665 - 663}, 

r2 = Co 1{62 — 63,62 -f 63 — 266, 62-63-1- 666 “ 669}, 

Fa = Co 1{63 — 64, 63 -f 64 — 267, 63 — 64 -|- 667 — 6610}, 

r4 = Co1{62 - 65,62 - 465 + 368}, r5 = Col{63 - 65, 64 - 67}, 

8 



re = col{e 3 - 466 + Seg, 64 - 4 e 7 + 3eio}, 

Uo{h) = Ee{n{hfp{n2 - ^^i)) + - n'^pn2, 

Hi = efQiei - e^Qie 2 + e|’( 52 e 2 - 64(5264, 


U 2 

ns 

He 


A^[hiRi + /112-R2 + 7(^1 ~ l)‘S'i + 7(^12 ~ 1 ) 5 ' 2 ]^, 
rfi?i(/ii)ri, n4 = 


2 (/'^, + i) rj5i(/ti)r4,n6 = 2rp2r5 + 4ri^52r6, 

til 


r2 

T 

R 2 X 


r2 

r3_ 


x^ R 2 


r3_ 


Riihi) = diag{Ri,3ciihi)Ri,5c2{hi)Ri}, 


R2 = diag{i?2, 3R2, 5R2}, 

Si{hi) = diag{5i,2c3(/ii)5i},52 = diag{52,52}, 


where Ci{hi) = I, i = 1,2,3, if /ii = 1 and ci{hi) = {hi + — 1), C 2 (/ii) = {hi + l)(/ii + 

2)y{{hi - l){hl + 11)), C3{hi) = {hi + 2)/{hi - 1) if hi > 1. 

The following reciprocally convex comhination inequality ifTTl will he used in the proof of our results. 


Lemma4: For given matrices Ri G S^, R 2 G matrix X G satisfying 

the inequality 


0 


1 - 

1 _ 


> 

0 J^R2_ 


* R 2 


Ri 

* 


X 

R2 


> 0 , 


holds for all a G (0,1). 

Proof: An elementary proof is derived from the fact that for any positive scalars o, b, the inequality 

—h -^ (\/o + ‘ s / h ' y ^ ^ a + 6 + 2c 

a 1 — a 

holds for all a G (0,1) and scalar c subject to ah > (?. ■ 

Theorem 1: Assume that there exist symmetric positive definite matrices P G Qi,Ri,Si G S+, 
i = 1, 2, and a matrix X G R3’^><3n following LMIs hold for h G {/ii, ^ 2 } 


i?2 ^ 

* R 2 


> 0 , 


2 6 

n(/i) = no(/i) + ^ Hi - ^ iij < 0 . 

i=l j=3 

Then system (l20l) is asymptotically stable for any time-varying delay h{k) G [hi,h 2 ]- 


( 21 ) 

( 22 ) 
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Proof: Consider the following LKF 

k-l 

=5F {k)Px{k) + ^ x'^(s)(5ix(s) 

s=k—h\ 

k—h\ — l —1 k—1 

+ ^ x^{s)Q2x{s) + hi Y, Y1 Ax^{i)RiAx{i) 

s=k—h 2 s=—hi i=k-\-s 

—hi — 1 k—1 

+ '•1. E E Ax'^ {i)R2Ax{i) 

s=—h 2 i=k-\-s 
—1 s k—1 

^ ^ ^ ^x^ij)SiAx{j) 

s=—hi i=—hi j=k+i 
—hi — 1 s k—1 

+ ^x'^ij)S2Ax{j), 

s=—h2 i=—h2 j=k+i 
k—1 k—hi — 1 —1 k—1 

where x(/i:) = col{x(fe), Y1 Y1 3;(s), Y1 Y1 a^(*)} and denotes the segment {x(A:) : 

s=k—hi s=k—h 2 s=—hi i=k+s 

k G Z[-/i2,0]}. 

The previous functional is positive definite due to the assumptions of Theorem [T] Now, we employ our 
newly derived inequalities in Lemma[2]and Lemma[3]in hounding Note at first that x{k + l) = 

(n(/i) — Qi)Co{k), x{k) = {Q{h) — Q, 2 )Co{k) and A{x^{k)Px{k)) = {x{k) + x{k + 1)1^PAx{k). Then 


we have 


AF(xI'=l) = Co^ {k) (no(/r) + Hi + U 2 ) Co{k) 

k—1 k—hi — 1 

— hi A^x{s)RiAx{s) — hi 2 A’^x{s)R 2 Ax{s) 


s=k—hi 
— 1 k+s 


s=k—h2 

—h\ — l fc+5 


(23) 


- E E A^x(z)5iAx(f) — E E A^x(z)52Ax(z). 

s=—hi i=k—hi s=—h2 i=k—h2 

Note that, the following equality 

b s b b b 

Y = ik-a + 2)Y v{s) YY 


(24) 


s=a i=s 


holds for any sequence v : Z[a, 5] M”. Using (|2^ in presenting ^3 defined in Lemma [2j we have 


k-l 


-hi Y A^x{s)RiAx{s) <-Co{k)TjRi{hi)TiCo{k). 

s=k—hi 

Similarly, the second summation term of (|2^ can he hounded hy (fT^ and Lemma |4] as follows 

k—hi — 1 


( 25 ) 


-hi2 Y^ A’^x{s)R2Ax{s) < 

s=k—h2 


f^l2 


h — h-\ 


C^{k)T^R2T2Coik) 
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hl2 


h2-h 
= -Co (k) 

^ “C(f (k) 


Co' (fc)r3i?2)r3Co(A:) 


1- 

(M 

1_ 

T r 

rg 


-1 

bO 

_1 

T r 

1- 

CO 

1_ 



hi2 

h—hi 

0 


R 2 0 



r2 


r3_ 


R 2 X 
* R 2 


Co{k). 


Co{k) 





r3_ 


Note that, when h = hi and h = h 2 then r 2 Co(^) = 0 and r 3 Co(A:) = 0, respectively, and thus the last 
inequality is still valid. Therefore 

k—hi — 1 

-hi2 X'^x{s)R2Axis)<-Cjik)U4Coik). (26) 

s=k—h2 

By Lemma [3] we have 

— 1 fc+s 

- Ax^{i)SiAx{i) < -C^{k)U^Co{k). (27) 

s=—hi i=k—hi 

Now, we employ ([T^ to hound the last term in (l23l) . To do this, note at first that 

—hi — 1 k-\-s —hi — 1 k-\-s 

E E Ax^{i)S2Ax{i) > E E Ax"’" {i)S2Ax{i) 


s=—h2 i=k—h2 


s=—h i=k—h 
—h—1 k-\-s 


+ E E Ax"’" { 1)82 Ax {i) 

s=—h2 i=k—h2 

Then, hy applying ([T^ and rearranging the obtained results we get 

—hi — 1 k-\-s 

- ^ ^ Ax'’"{i)S 2 Ax{i) < -Coik)UeCoik). 

s=—h2 i=k—h2 

It follows from (I23l)-(l2^ that 

AL(x["l)<Co^(fe)n(/r)Co(fe). 


(28) 


(29) 


The matrix n(/i) is an affine funcfion h, and thus, n(/i) < 0 for all h £ [hi, /i 2 ] if and only if n(/ii) < 0 
and n(/i 2 ) < 0. Therefore, if (l22b holds for h = hi and /i = /12 then, from (l29l) . A17(x[^]) is negative 
definite which ensures the asymptotic stability of system (l20l) . The proof is completed. ■ 

Remark 6: In Theorem [H a full 3n x 3re matrix X is used in the reciprocally inequality to improve 
the upper bound of delay. However, to reduce the number of decision variables, we can use the matrix 
X of the form X = diag{Ai, X 2 , A 3 }, where A* E f = 1, 2,3. 
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B. Examples 


Example 1: Consider system (l20l ) with the matrices taken from the literature 



0.8 

0.0 


-0.1 

0.0 

A = 

0.05 

0.9 

) — 

-0.2 

-0.1 


The obtained results and comparison to most recent results in the literature are given in Table I and 
Table II. It is worth noting that, thanks to our new summation inequalities proposed in Lemma |2] and 
Lemma[3l Theorem [T] clearly delivers significantly better results than the existing methods in the literature. 
Especially, our method requires less decision variables than the proposed conditions in ll22ll . lf24ll . ll25l 
while leading to much better results. 


TABLE I 

Upper bounds of /12 for various hi in Example 1 


hi 

2 

4 

6 

10 

15 

20 

25 

30 

NoDv 

(Proposition 1) 

17 

17 

18 

20 

23 

27 

31 

35 

8 n^ + 3n 

nn (Theorem 1) 

18 

18 

19 

20 

23 

26 

30 

35 

3.5n^ + 3.5n 

(271 (Theorem 3.1, ( = 4) 

20 

21 

21 

22 

24 

27 

29 

34 

9.5n^ + 5.5n 

(721 (Theorem 4, ( = 4) 

21 

21 

21 

22 

24 

27 

31 

35 

9.5n^ + 5.5n 

1221 (Theorem 2) 

22 

22 

22 

23 

25 

28 

32 

36 

27n^ + 9n 

1241 (Theorem 2) 

22 

22 

22 

23 

25 

28 

32 

36 

23n^ + 7n 

(251 (Theorem 1) 

22 

22 

22 

23 

26 

29 

32 

36 

19n'^ + 5n 

Remark [6] 

24 

26 

27 

30 

32 

33 

35 

39 

14n^ + 5n 

Theorem [T] 

26 

27 

28 

31 

34 

35 

36 

39 

20n^ + 5n 


NoDv: Number of Decision variable 


TABLE II 

Upper bounds of /12 for various hi in Example 1 

hi 

I 

3 

5 

7 

11 

13 

NoDv 

(H 

12 

13 

14 

15 

17 

19 

9iA + 3n 

GH 

17 

17 

17 

18 

20 

22 

13n^ + 5n 

dD 

17 

18 

19 

21 

25 

25 

Dv« 

(261 

20 

21 

21 

22 

23 

24 

10.5n^ + 3.5n 

Thm.[T] 

26 

, .'2 , /I 

27 

28 

29 

32 

33 

14n^ + 5n 


Dv* = (/12 + lfn^/2 + {h 2 + 2)n/2 
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Example!: Let us now consider a practical satellite control system Il28ll . The dynamic equations are 
as follows 


JMt) + - ut)) + k{ei{t) - 92{t)) = uit), 

J2Ht) + - Ht)) + kiei{t) - d2{t)) = 0, 


(30) 


where Jj, i = 1,2, are the moments of inertia of the two bodies, / is a viscous damping, A; is a torque 
constant, 9i{t) are the yaw angles for the two bodies and u{t) is a control input. The following parameters 
are borrowed from ll28l : Ji = J 2 = 1, A: = 0.09,/ = 0.04. Let Xi{t) = 9i{t), Xi+ 2 {t) = 9i{t), i = 1,2. 
By choosing a sampling time T = 10 ms, system (l30l ) can be transformed to the following discrete-time 
system Il2^ 


x{k + 1) = Ax{k) + Bu{k) (31) 

where 


1 

0 

0.01 

0 


0 

0 

1 

0 

0.01 

,B = 

0 

-0.009 

0.009 

0.9996 

0.0004 

0.01 

0.009 

-0.009 

0.0004 

0.9996 


0 


A delayed state feedback controller is designed in the form u{k) = Kx{k — h{k)), where h{k) is 
time-varying delay belonging to the interval [hi, / 12 ]. For hi = 1, it was found that with the controller 


gain K = 


0.1284 -0.1380 -0.3049 0.0522 


, Theorems 1 and 2 in 11221 give the upper bounds of 
/i 2 as 129 and 135, respectively, which are larger than 98 delivered by the results of lITSl . We apply 
Theorem [T] for Ad = BK, it is found that the closed-loop system remains asymptotically stable for the 
time-varying delay h{k) G [1,170] which shows a clear reduction of the conservatism. To demonstrate 
the effectiveness of the obtained result, a simulation with h{k) = 1-|-169| sin(A:7r/2)| and initial condition 
[2—1 0.2 — 0.5]^ is presented in Fig. 1. It can be seen that the state trajectory converges to zero as 
shown by our theoretical result. 


V. Conclusion 

In this paper, new summation inequalities in single and double form have been proposed. By employing 
the newly derived inequalities, improved stability conditions have been derived for a class of discrete-time 
systems with time-varying delay. Provided examples and comparisons to the most recent results found 
in the literature show the potential and a large improvement on the stability conditions deliver by the 
approach proposed in this paper. 
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10000 


15000 


Fig. 1. Responses of the satellite system with h{k) = 1 + 169| sin(fc7r/2)| 


References 

[1] R. Sipahi, S.-I. Niculescu, C.T. Abdallah, W. Michiels, and K. Gu. Stability and stabilization of systems with time delay. 
IEEE Control Systems, 31(1): 38-65, 2011. 

[2] E. Fridman. Introduction to Time-Delay Systems: Analysis and Control, Birkhauser, 2014. 

[3] E. Fridman, U. Shaked, and K. Liu. New conditions for delay-derivative-dependent stability. Int. J. Control, 45(11): 2723- 
2727, 2009. 

[4] L.V. Hien, and V.N. Phat. Exponential stability and stabilization of a class of uncertain linear time-delay systems. J. 
Franklin Inst., 46(6): 611-625, 2009. 

[5] E. Fridman. A refined input delay approach to sampled-data control. Automatica, 46(2): 421-427, 2010. 

[6] A. Seuret. A novel stability analysis of linear systems under asynchronous samplings. Automatica, 48(1): 177-182, 2012. 

[7] A. Seuret, and F. Gouaisbaut. Wirtinger-based integral inequality: Application to time-delay systems. Automatica, 49(9): 
2860-2866, 2013. 

[8] L.V. Flien, and H.M. Trinh. A new approach to state bounding for linear time-varying systems with delay and bounded 
disturbances. Automatica, 50(6): 1735-1738, 2014. 

[9] Y. Li, K. Gu, J. Zhou, and S. Xu. Estimating stable delay intervals with a discretized Lyapunov-Krasovskii functional 
formulation. Automatica, 50(6): 1691-1697, 2014. 

[10] P.G. Park, W.I. Lee, and S.Y. Lee. Auxiliary function-based integral inequalities for quadratic functions and their applications 
to time-delay systems. J. Frankl. Inst., 352(4): 1378-1396, 2015. 

[11] P.G. Park, J.W. Ko, and C. Jeong. Reciprocally convex approach to stability of systems with time-varying delays. 
Automatica, 47(1): 235-238, 2011. 


14 

















[12] Z. Feng, J. Lam, and G.Fl. Yang. Optimal partitioning method for stability analysis of continuous/discrete delay systems. 
Int. J. Robust Nonlinear Control, 25(4): 559-574, 2015. 

[13] C. Peng. Improved delay-dependent stabilisation criteria for discrete systems with a new finite sum inequality. lET Control 
Theory Appl, 6(3): 448^53, 2012. 

[14] A. Seuret, and F. Gouaisbaut. Complete quadratic Lyapunov functionals using Bessel-Legendre inequality. European Control 
Conference (ECC’14), Strasbourg, France, 2014. 

[15] L.V. Hien. A novel approach to exponential stability of nonlinear non-autonomous difference equations with variable delays. 
Appl. Math. Lett., 38: 7-13, 2014. 

[16] E. Fridman, and U. Shaked. Stability and guaranteed cost control of uncertain discrete delay systems. Int. J. Control, 
78(4): 235-246, 2005. 

[17] Y. Fie, Q. G. Wang, L. Xie, and C. Lin. Further improvement of free-weighting matrices technique for systems with 
time-varying delay. IEEE Trans. Autom. Control, 52(2): 293-299, 2007. 

[18] B. Zhang, S. Xu, and Y. Zou. Improved stability criterion and its application in delayed controller design for discrete-time 
systems. Automatica, 44(11): 2963-2967, 2008. 

[19] X.Y Meng, J. Lam, B.Z. Du, and H.J. Gao. A delay-partitioning approach to the stability analysis of discrete-time systems. 
Automatica, 46(3): 610-614, 2010. 

[20] H. Shao, and Q.L. Han. New stability criteria for linear discrete-time systems with interval-like time-varying delays. IEEE 
Trans. Autom. Control, 56(3): 619-625, 2011. 

[21] C.-Y. Kao. On stability of discrete-time LTl systems with varying time delays. IEEE Trans. Automatic Control, 57(5): 
1243-1248, 2012. 

[22] O.M. Kwon, M.J. Park, J.H. Park, S.M. Lee and E.J. Cha. Stability and stabilization for discrete-time systems with time- 
varying delays via augmented Lyapunov-Krasovskii functional. J. Frankl. Inst., 350(3): 521-540, 2013. 

[23] L.V. Hien, N.T. An, and H. Trinh. New results on state bounding for discrete-time systems with interval time-varying delay 
and bounded disturbance inputs. lET Control Theory Appl, 8(14): 1405-1414, 2014. 

[24] S.H. Kim. Further results on stability analysis of discrete-time systems with time-varying delays via the use of novel 
convex combination. Appl. Math. Comput., 261: 104-113, 2015. 

[25] P.T. Nam, P.N. Pubudu, and H. Trinh. Discrete Wirtinger-based inequality and its application. J. Erankl. Inst., 352: 1893- 
1905, 2015. 

[26] A. Seuret, F. Gouaisbaut, and E. Eridman. Stability of discrete-time systems with time-varying delays via a novel summation 
inequality. IEEE Trans. Autom. Control, doi: 10.1109/TAC.2015.2398885. 

[27] J. Zhang, C. Peng, and M. Zheng. Improved results for linear discrete-time systems with an interval time-varying input 
delay. Int. J. Syst. Sci., (2014) doi: 10.1080/00207721.2014.891674. 

[28] H. Gao, X.Meng, T. Chen, and J. Lam. Stabilization of networked control systems via dynamic output-feedback controllers. 
SIAM J. Control Optim., 48(5): 3643-3658, 2010. 


15 



